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ABSTRACT 

We show that the large-scale cosmic ray anisotropy at ~ 10 TeV can be explained by a modified Compton- 
Getting effect in the magnetized flow field of old supernova remnants. This approach suggests an optimum 
energy scale for detecting the anisotropy. Two key assumptions are that propagation is based on turbulence 
following a Kolmogorov law and that cosmic ray interactions are dominated by transport through stellar winds 
of the exploding stars. A prediction is that the amplitude is smaller at lower energies due to incomplete sampling 
of the velocity field and also smaller at larger energies due to smearing. 



1. INTRODUCTION 

In the past few years, a number of cosmic ray (CR) experi- 
ments have convincingly shown that the arrival directions of 
Galactic CR particles are not fully isotropic in the sky: A 
large-scale anisotropy with a dipole and quadrupole compo- 
nent was first observ ed in the northern hemisphere at an en- 
ergy of several TeV ( Amenomori et al. 2006; Guillia n et al. 
2007) |Abdo et al.|2008[|2009| ). The same structure was later 
found in the southern hemisphere at a mean energy of 20 TeV 
( Abbasi et al. 2010| >. The large-scale anisotropy /, observed at 
several to tens of TeV energies, is of a level/ - 10" 4 to 10" 3 . 
The results are summarized in Table [TJ 
Apart from the large-scale anisotropy, several experiments 



TABLE 1 

The ob served large-scale CR a nisotropies. Referen ces are: 

[1] -IAbbasi et al.1(|2010|i ; [2 ] -|Abbasi et al.M2012> ; [3] - 
IAmenomori et al.|i|2U06|i; [4] -|AMi)0 et al.H2008| i; [5] - |Abdo | 



degree scale and with an intensity of the order of < 10 -4 ( Ab 
basi et al.|201 1 1. While the deficit has a significance of 6.3 a, 
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showed that at the same energy scale smaller angular scale 
excesses and deficits do exist, with extensions from a few de- 



grees in the sky up to about 20 degrees in one direction ( Abdo 
|et al.|2 009 ; Abbasi et al. 2011). Most recently, it was shown 
that the large-scale anisotropy which is present at TeV ener- 
gies vanishes in the southern hemisphere at around 400 TeV 
and a new component emerges instead, a clear deficit of a 20 
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a clear excess is not yet distinguishable from current statistics 
(Abbasi et al. 2012). What is important about this result is 



that it is clearly not a simple dipole component, which would 
be just symmetric in excess and deficit, but the anisotropies 
clearly display disjunct excesses and deficits. 
First we may have to ask, why we should expect isot ropy of 
cosmic rays at all: the essential answer was given by Schluter 
|and Biermann | ( fl950| ): Magnetic fields get strengthened until 
they scatter cosmic rays into n ear prefect isotropy (also found 
by Hanasz et al.||2004" 2009 1. On this basis they estimated 
the strength of magnetic fields to be of order 4 /nG, an es- 
timate, which has held up remarkably well (e.g. Beck et al. 
1996). Thus we can expect isotropy in the reference frame 
of the local interstellar me dium, which is ve ry well coupled 
to the magnetic fields (e.g. |Appenzeller][T974|, a connection 
which is kept by instabilities (Parker 1966 ). However, the Sun 
was only coupled to the interstellar medium at birth, 4.5 bil- 
lion years ago: Interaction with massive interstellar gas clouds 
slowly increases the pecu liar velocity of stars with age (e.g. 
|Julian|1967j |Wielen|1975| ), and so at the age of stars like the 
Sun the increase in peculiar velocity is expected to be about 
40 km/s, which is ac tually somewh at more than deduced from 
3D-observations ( |Reid et al.|20 09 ). Therefore, within such a 
velocity the cosmic ray anisotropy should be small relative to 
the average interstellar medium around the Solar system. 
In this paper, we discuss the possible origin of the large-scale 
structure anisotropies at — 10 TeV with some comments on 
larger and smaller energies. 

2. COSMIC RAYS AND THEIR GALACTIC 
PROPAGATION 

Cosmic ray particles are believed to be injected into the ISM 
by supernova explosions, as soon as their shocks slow down 
sufficiently to release the population of energetic particles; 
this may happen for progenitor stars of modest mass as explo- 
sions directly into the interstellar medium, or for very mas- 
sive stars as explosions into their stellar winds. The distinc- 
tion between these two types o f supernova explosions (Stanev 
2012| and earlier papers) allows to 



et al. 19931 Nath et al. 



interpret th e cosmic ra y positron and electron enhancements 
( |Biermann~e t al. 2009 ), the WMAP haze of high frequency ra- 
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dio emission as well as the 511 keV annihilation line near the 



-3.26±0.06 



Galactic Center (Biermann et al. 2010a i, as well as th e upturn 
in the cosmic ray spectra of nuclei (Bierma nn et al.| [2010b). 
Finally, it allows to understand the KASCADE-Grande spec- 
tra of cosmic ray par ticles at energies beyond 10 15 eV (Bier- 
mann and de Sou"za||2012[ l, as well as allowing to use the 
cosmic ray particl es as injection seeds for the ultra-high en- 
ergy cosmic ray s (|Gopal-K rishna et al.||2010 Biermann and 
de Souza 2012). Then the cosmic ray particles meander, af- 



up to TeV energies ( Wiebel-Sooth and Biermann 



ter injection, through the interstellar medium (ISM) subject 
to the irregularities of the magnetic field. In a general flow 
not too far from equipartition, the spectr um of the irregular- 
ities should be roughly Kolmogorov (e.g. Kolmogoro v| 1 94 1 1 
Sagdeev 1979 Goldstein et al.|1995 1. This has been shown to 



be consistent with observations (Sprangler and Gwinn 1990| >. 
This has been demonstra ted to oe consis t ent with the low 
level of anisotropies (e.g. IBiermannll 19931 IBlasi and Amatol 
|2012a|b| ). 

2.1. Interaction test with cosmic rays 

However, it is well known, that at the same time, as parti- 
cles scatter around they also interact, and so the spectrum of 
the spallation secondaries, like Boron, compared with Car- 
bon, should also reflect the same energy depe ndence: This i s 
not correct, since the B/C ratio scales as E~° 54 (Ptuskin 1999); 



this is inconsistent with a Kolmogorov spectrum. Ho wever, it 



could be c onsistent with a Kraichnan spectrum k~ 3 ^ 2 ( Kraich 
|nan|| 196 51): However, Kraichnan turbulence implies a lower 
dimensionality due to a dominant ma gnetic field, wh ich is 
not given in the interstellar medium ( |Beck et al.|1996| l. The 
contradiction can be resolved: The massive stars, that pro- 
duce most of the CR Carbon before they explode, are Wolf 
Rayet stars. These stars have a powerful wind and eject most 
of their zero-a ge-main- sequence mass into the wind, before 
they explode (Prantzos 1984). This wind forces most of the 
stellar mass as well as a large amount of ISM material into 
a thick shell. Therefore at the time of the explosion, typi- 
cally more than half of the original stellar mass is contained 
in the wind and its shocked shell, composed of old stellar and 
ISM material. There, th e CRs themse l ves excite a spe ctrum 
of magnetic turbulence ( |Bell|1978a|b[ ). [Biermann ( 1998) cal- 
culated the energy dependence of the escape yielding for the 
B/C ratio an energy dependence of E~ 5 / 9 based on the e xcita- 
tion of th e irregularities by the cosmic rays themselves (Bell 
|1978a|b|) is fully con sistent with the observed spectral depen- 
dence ( |Ptuskin|1 999 ). This suggests that most of the CR inter- 
action happens in these shells. A recent example of observing 
thi s cosmic ray int eraction directly in 7-ray data was shown 
by Berezhko et al. (2004), who argued directly for interaction 
in a wind-environment. It also implies that there is a minimum 
interaction, or in other words, a finite minimum path-length of 
interaction, stemming from the fact that the particles at very 
high energies do not scatter through the wind-shell, but con- 



vect ( |Biermann|1993l|Biermann et al.|200~T]|Nafh et al.|2012| i. 
In this latter point of a finite residual path length the model is 
consistent w ith the results by the Tracer experiment (Ober- 
|meier|[20TT| ), but needs larger statistics data for final confir- 
mation. There is another test for this argument from a com- 
parison of the electron and proton spectra, which should have 
the same injection slope at TeV energies if arising from the 
same shock mechanism in the same source class: At ener- 
gies well above 10 GeV the CR electrons have a spectrum of 



1999) and clearly have been steepened by synchrotron and 
inverse Compton losses ( Kardashe vj 1 962| l, and so their injec 
tion spectrum is £~ 2 26=l=Q ut '. Comparing this with the CR pro 



tons (CREAM: Yoon et aT]|201 1 1, which give a spectrum of 
about E~ 1M±Qm near TeV energies, the difference to the cor- 
rected CR electron spectrum gives the energy dependence of 
the diffusive escape, £ , -°- 40 ± 06 5 fully consistent with a Kol- 
mogorov law. A spectrum of close to E~ 8 / 3 matches the p re- 
diction for wind-SN CRs (see also Biermann et al. (2010b) on 
the match to the new CREAM data ( [Yoon et al.||201l| X On 
the other ha nd, the arguments for t he positron fraction rising 
with energy (Bierma nn et al.|2 009 ) support the point of view 
that wind-SNe are more important for observed CR-electrons 
at these en ergies; the pred icted ratio of CR-positrons to CR- 
electrons pJiermann et al. 2009)) has recently been confirmed 



again by|Ackermann et aL 



(2012). Therefore, in the following 



we will use the concept, that for interstellar medium propaga- 
tion the relevant spectrum of irregularities has a Kolmogorov- 
type spectrum. 

2.2. Isotropy and anisotropy 

However, before we move on to discuss the anisotropies we 
should note, that the observation of a near-isotropy is al- 
most more interesting: The isotropy has been used already 
early t o argue that magnetic fields are required to isot ropize 
them (|L. Biermann and Schluter|195T]|L. Biermann|1953a|b| 
|L. Biermann and Davis|1958||1960||Ginzburg and S yrovatskii 

1964), and so their strength has been estimated remarkably 
close to what has been measured since (e.g. Beck et al. 1996). 
But what had never been clear is whether the cosmic ray refer- 
ence frame is really our frame, the Solar system frame, or any 
specifically defined local frame, with reference to local stars 
for instance. In order to discuss the origin of anisotropies we 
need to consider what we actually observe: Cosmic ray parti- 
cles scatter with a mean free path dependent on the spectrum 
of magnetic irregularities; so in any direction we observe the 
surface of last effective scattering, and it is critically impor- 
tant to note, that this surfaces itself has a depth equivalent to 
its distance, so for any specific particle the distance is defined 
by the length scale over which it forgets its original direction, 
and comparing different particles that distance is a broad dis- 
tribution with a width similar to its distance, as indic ated by 
the exponential pa th-length distribution description ( |Garcia-| 
Muno z~et al.|1987] ). It immediately follows that looking into 
different directions, by about a radian, these patches of the 
surface of last effective scattering may well be uncorrected, 
with of order 4 patches possible, one of which will be domi- 
nant, one will be second, and the rest will produce what looks 
like noise, unless our statistics become incredibly good. Of 
course if one recent SN-explosion completely wipes out all 
traces of previous other SN-explosions then we will be able 
to detect that in many other traces of recent SN-explosions, 
possibly even with light-echoes; so individual sources could 
produce clearly excesses of smaller angular extent. Similarly, 
magnetic enhancements in specific features, like the heliotail, 
could produce deficits of smaller angular extent. 

Figure 1 shows a schematic picture of the model presented 
here, emphasizing the role of the surface of last effective scat- 
tering 



Using a diffusive des cription of cosmic ray propagation Blasi 
|and Am ato (2012ab) have modelled the possible cosmic ray 



Cosmic ray transport and anisotropics 



3 



Surface of last effective scattering 
(dependent on energy) 




FIG. 1 . — Schematic geometry of the surface of effective last scattering. We 
note that this surface is itself quite "thick" with a range of scattering distances 
of the same scale as the scattering path itself 

anisotropics, and have found that the anisotropy could be 
quite large. 

Here we propose to take a slightly different route, comparing 
the influence of both diffusion and velocity fields, always em- 
phasizing that the surface of last effective scattering is "thick", 
as noted; this entails that the higher the energy, the larger the 
scale, but also the larger the smearing over all lower scale 
spatial variations. We will focus on the large angle anisotropy 
near ~ 10 TeV. 

The escape from the Galactic disk can be described as a ran- 
dom walk with the step size a function of energy ( |Chan-| 
|drasek har 1943 ). Another related description is with the ex - 
ponential path-length distribution ( |Garcia-Munoz et aL[ 1987 ). 
The escape time from the thick Galactic disk (| Beuermann 



|et al.|1985||Ferrando|1993]|Brunetti and Codino|2000) of an 
approximate half thickness of d = 1 -2 kpc is of order 10 7 yrs, 
and using the argument above depends on energy as E' 1 ^ . 
Considering a possible heavy composition of Galactic CRs at 
1Q18.5 e y we w jjj sca j e iq proton energies at 10 17 eV. Under 
the cautious assumption that the thickness H is of order three 
times the mean free path at the maximal energy we obtain a 
mean free path A m f p of order: 



Amfp = -r H 



10 17 eV 



1/3 



10 1 



VGeV j 



1/3 



(1) 



Since at that level in vertical direction, about 1-2 kpc, we 
have a tran sition to a Galactic magnetic wind ( |Everett et al. 
|2008| |2010[ |Everett and Z weibel 20TTJ, which also scatters 
particles quite effectively due to the larger scale and the 1 jr- 



behavior of the magnetic field (Parker 1958 ; Biermann and de 



Souza[ 
The scale 



120T21 , 
:ale of 



this ought to be enough to ensure near isotropy. 



order 20 pc is then reached for particle energies 
of about 10 TeV with an energy dependence as £ ,1//3 , and a 
scale of order 70 pc is reached at about 400 TeV. 
So, let us first consider the various sources of anisotropics, 
first looking at source distributions, and second at velocity 
fields. 

2.2.1. Sources 



The first effect is an overa ll gradient of cosmic ray particle s 
in a disk galaxy like ours (Ginzbu rg and Syrovatskii ||1964 l; 
however, a disk galaxy is usually a spiral, with some variation 



across the spiral arms (e.g. Beck and Hoern esfl996 l, so that 
the key radial gradient has to be considered along the spiral 
arms, so effectively increasing the radial scale of comparison 
by {cos 9}~ l . where 8 is the angle of the spiral arms with re- 
spect to a circle. This implies that we expect anisotropy of 
order one permille, by dividing 20 pc by several times 8 kpc, 
but symmetric, which is not seen in the data. This probably 
implies that this effect is washed out out by other influences. 
Detailed observations a t several wavelength in radio and in- 
frared clearly show (e.g. Tabatabaei et al. 2007), that the local 
variations on scales of a few tens of pc and larger dominate the 
unevenness of the non-thermal radio-emission, and so proba- 
bly also of the cosmic ray distribution. At much higher ener- 



gies th is effect may have been detected (see e.g. Teshima et 
al.|2001 1, but even there, source regions may dominate what 



is seen. 

The second effect is due to the diffusion stra ight out from the 
disk into the halo Galactic magnetic wind ([Everett et al.|2008~ 



|20101 [Everett and Zw eibel 2011). However, this effect is of 
second order due to the symmetry, and depends critically on 
where the mid-plane is for this flow and diffusion, which we 
do not know. Assuming for lack of a good number for the dis- 
tance from Earth to mid-plane of the flow-field about 40 pc, 
and for the vertical scale of the flow again 1 .5 kpc, suggests an 
asymmetry of order 6 x 10~ 4 , no longer visible at any length 
even approaching 40 pc. This again is clearly not dominant in 
what we observe. 

Finally, the obvious most recent sourc e, an idea which 
has been explored very many times (e.g. |Volk et aL[|1988| 
Teshima et aT[[200l] |Erlykin and W olfendale[ [2006[ lYuksel] 
et al.|2009| ). One difficulty with such an idea is that there is no 



clear spectral signature of any recent source of cosmic rays; 
on the other hand, such a source must exist, just statistically. 
The question is, what is its signature? The supernova rate in 
our Galaxy has been estimated to be of order 1 per 30 years, 
with a n uncertainty towards longer times of a factor of 3 (see 
e.g. Bierm ann and Cassinelli|1 993 ; Biermann et a l.|1995| ); av- 
eraging this uncertainty blithely suggests that one supernova 
in 50 years, to within a factor of 2, might be a good initial 
estimate. The probability to then have a supernova explode 
within a certain radial distance over a certain time, at our dis- 
tance from the Galactic Center, can then be estimated. Since 
the supernova rate is reduced at our distance from the Galactic 
Center with respect to the average can be crudely estimated to 
be of order 3, so that in our part of the Galaxy, the supernova 
rate may be about one in 150 years. It follows that within a 
distance of r = Ri .3 20 pc from us the time scale between su- 
pernovae is of order 10 7,6 /?7 2 yrs. 

This, then, allows to derive the diffusion coefficient derived 
above as 10 29 8 cm 2 s~ 1 at 10 TeV; this is slightly larger than 
the estimate given by |Blasi and Am ato (2012a b ), who use a 
different concept of cosmic ray propagation, but our number 
is close to older estimates. The uncertainties are dominated by 
systematics, since different lines of reasoning have been used 
to derive the coefficient as well as its energy dependence. 
This implies that over a distance of 20 pc the diffusive time- 
scale is of order 10 2,4 ,f? 2 3 yrs. Therefore the dispersion of 
any original cloud of CR-particles is extremely effective, with 
the typical time scale between supernovae taken from above 
we obtain a radial range of many kpc numerically, and a di- 
lution relative to the average cosmic ray energy density of 
about 10~ 5 4 , so the typical anisotropy expected from super- 
novae should be quite small. This predicted anisotropy scales 
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with the reference distance as /?j 3 3 , so for a mean free path at 
10 TeV larger by a factor of two, well within the uncertain- 
ties, we obtain a predicted cosmic ray anisotropy of 10" 45 , 
getting close to numbers observed. This clearly implies that a 
recent supernova explosion cannot be immediately discounted 
as a contributin g origin for cosmic r ay anisotrop ies, confirm- 
ing (e.g. [Vblk et al|l98"8) |Teshima et al.pOOTj |Erlykin and 
|Wolfendale|2006t|Yuksel et al.|2(W)| and others). 

2.2.2. Residual flow fields 

However, when we consider the flow field produced by super- 
nova explosions, then the numbers change, since the flow is 
just mixed slowly after decaying to subsonic and sub-Alfvenic 
velocities. 

The simplest second reason for cosmic ray particles to show 
anisotropy is a general movement of our fram e of reference 
with respect t o that of the cosmic ray system (Compton and 
|Getting|1935| l. 

Data and magneto-hydrodynamic simulations suggest that the 
ISM is mainly driven by supernova explosions, phase transi- 
tions, cloud and star formation, radiation from young stars, 
instabilities and energy transport by magnetic fields and CR 
particles, and shear and outflow from the Galactic disk ( Cox 



1972 



et al. 



Parker 1966; McKe e and Ostriker| 19771 IBreitschwerdt 
19911)Lee et al.|20031|Hanasz et al.|2004] |Everett et al. 
2008| |2010[ |Everett and Zweibel||201 1| >. The hottest normal 
phase of the ISM is at about a de nsity of 3 x 10~ 3 cm" 3 and a 



temperature of a bout 4 x 10 6 K ( Snowden et al.|l997 Hagi 



hara et al.|201 l| l, confirming earlier expectations (Lagage and 
Cesarsky [T983). After a supernova explodes into the ISM, 



the expansion of the shock racing through the ISM is fast at 
first, accumulati ng evermore m aterial from the ISM, and then 
slows down ( |Sedov|1958[ |Cox|1972) , until the flow becomes 
subsonic and sub-Alfvenic; thereupon the flow coasts along, 
and basically ge ts slowly disorganiz ed by mixing and encoun- 
tering clouds. |Gaensler~eT al. (2011) confirm that the typical 
sonic Mach-number is low, of order 2 or less. At the density 
and temperature of this most tenuous phase the signal speeds 
will thus be around 160 km/s, and so the velocity fields of the 
late evolution of supernova remnants (SNRs) will run at or a 
bit below these velocities, but later decay rather slowly. The 
average density across all m edia in the 200 pc thick layer of 
the ISM is of order 1 cm" 3 (Cox||2005|). Using this average 



density we obtain a velocity scale of order 100 km/s at a scale 
of about 30 pc, using the simple expressions of |Cox| ( fl972| l. 
The associated time scale is of order 2 x 10 5 yrs. Thus, the 
SNR is no longer supersonic or super-Alfvenic in the hot 
medium. Observational data on pulsar acti vity and their as - 
sociated SNRs support the scale of 30 pc ( |Braun et al.j 1989 >, 
and a time scale for the powering of 2 x 10 4 yrs, less than 
2 x 10 5 yrs, so consistent with the numbers suggested here for 
the hot interstellar gas. Similarly, the confinement time for 
CR particles in decaying SNRs has been estimated also to be 
relatively short, of order 10 4 yrs (Berez hko and Volk 2004| l, in 
agreement with the arguments by Braun et al. ( 1989). Cosmic 
rays will travel most easily in that phase of the ISM, where 
the Alfven velocity is the highest, due to limiting the stream- 
ing instability at the Alfven velocity. Therefore, we might 
expect the CR particles that we observe to come most effec- 
tively through this phase. However, as the observations show 
( | Appenzeller|| 1 974} the magnetic fields permeate the clouds, 
and so it is more suitable to use the average density of 1 cm" 3 . 



So the old velocity field of SNRs having gone subsonic and 
sub-Alfvenic in the highest temperature medium corresponds 
to velocities of about 100 km/s and to length scales of about 
30 pc. 

Considering the scattering of CRs arriving at Earth, at such 
a distance from us when they first point in our general di- 
rection, that surface can be called the surface of last effec- 
tive scattering, one mean free path away. If that region has a 
general flow field, it will imprint an anisotropy upon the CRs 
coming to Earth. The observed CR anisotropy corresponds 
maximally to about 10~ 3 ; comparing with the anisotropy of 
3 x 10" 4 attributed to the motion of the Earth around the Sun 
of 30 km/s, this implies a velocity amplitude of about 100 
km/s. This matches the velocities of old SNRs after they be- 
come sub-s onic and sub-Alfvenic, u sing the Compton-Getting 
effect (Compton and Getting 1935), which connects velocity 
field v, CR spectrum Zi^'dE and CR anisotropy A/// av : 



— = -(p + 2)cos( 

hy C 



(2) 



It has to be noted, that this effect does not require that we sit in 
the center of a SNR: The effect is there as soon as we have co- 
herent motion at the surface of last scattering over about one 
radian in one direction laterally, and in depth by about the dis- 
tance itself, so Ar/r ~ 1; consistent with this the data suggest 
an effect which varies across the sky. Massive stars that later 
produce supernovae f orm in the cold disk which is onl y of or- 
der 200 pc full width ( |Cox and Smith|1974||Cox|2005) . Using 
the time scale from above of about 200,000 yrs yields a length 
scale about a factor of 2 larger, with the energetics adopted by 



Cox ( 1972). The irregular flow field will be dominated by a 



single most recent supernova up to this scale combined with 
several supernovae that are somewhat older. However, us- 
ing this time scale we need to go back to our earlier question 
and ask, what the source anisotropy could be at such a time 
scale, and the answe r is 10" 3 6 , so in agreement with Blasi 
|and Ama"to| ( |2012a[bl ) not completely negligible compared to 
the influence of the flow itself. And also, just like the flow 
the directionality is not isotropic, the only key difference is, 
that the source contribution can only be positive, but the modi- 
fied Compton-Getting effect can be either positive or negative; 
the observations suggest at 400 TeV that the major effect is a 
deficit, consistent with a modified Compton-Getting effect. 
To summarize in yet another way, for any scale L <H, where 
H is the cosmic ray scale-height (Bier mann et al.pOOl 
have the inequality 



— < — 

^signal 



(3) 



where V s i g nai is either the speed of sound c s , the Alfven veloc- 
ity Va, or some combination thereof (the fast magneto-sonic 
speed). On the left side this is the cosmic ray diffusion time 
scale over a length L, and on the right hand side this is the 
convection time scale over the same length scale. At L = H 
the two time scales become equal, and convection takes over, 
and the loss from the disk is then driving a convective Galactic 
wind (B reitsch werdt et al.|1991[|Eve"rett et al.|2008"l|2010[|Ev-| 
|erett and"Zw eibel 2011 1 Bier mann et al.|2010a| l. However, it 
then also follows that convective flow is slower for all L < H, 
so lives longer, and can dominate cosmic ray scattering, and 
so is more viable to let us understand the observed anisotropy. 
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2.2.3. Magnetic field isotropy? 

There are several other effects which need to be considered in 
the propagation of CRs. First of all, the spectrum of irregular- 
ities may be well approximated by an isotropic Kolmogorov 
spectrum, but since the obs erved magneti c field distribution is 
never completely irregular (B eck etal.| 199 6) there can always 
be effects from this underlying anisotropy. Furthermore, even 
the irregularities themselves are never perfectly isotropic ei- 



ther (e.g.JMalkov et al.|20T0 
Lazarian 201T|i 



Lazarian and Desiati|2010 De- 



siati and Lazarian 201 1), and so such an assumption has to 
be taken with great caution, even though it may appear that 
it works relatively well. Finally, and perhaps most impor- 
tantly, the history of recent supernova explosions in the Solar 
neighborhood will give an imprint of irregularity just from the 
source distribution; however, since the time to escape is much 
larger than the time to replenish the cosmic ray population this 
should perhaps not be dominant until one gets to really high 

energies. 

The s ky maps in radio r otation measure (|O pperman n et al.] 
201 in radio emission (Berkhuij sen et al.||1971 1 and other 



wavelengths usually integrate over much larger distances, that 
it is difficult to be sure of a correlation with the detected 
anisotropics. There may be a correlation of s ome radio fea- 
tures like Spur 185- (Berkhuijsen et al.||197T} with features 
in Tibet, Milagro, SuperK and IceCube data; there may be a 
correlation with Loop I in IceCube data, and Loop IV in Tibet 
data; and finally, there may be a correlation with Spur 195+ 
in Milagro data. The deficit in the IceCube flux at 20 TeV is 
remarkably close to the center of Loop I. It needs however a 
careful analysis of the cosmic ray data and a detailed model- 
ing of the theory in order to investigate the causality of the 
signatures; today's radio data are much better (e.g. |Opper- 



|mann et al.|[20TT] |van Eck et aLlpOTT] |Pshirkov||201 1) , but 
what we need is yet more significant cosmic ray data. 
After subtracting the large angle fe atures it is possible to iden- 
tify small angular scale features (Abdo et al. 2008), w hich 
may related to Solar wind effects such as the heliotail (|Na- 



gashima et al. 1998 1 Drury and Aharonian 2008 Lazarian and 
Desiati 2010; Desiati and Lazarian 2011k or nearby stars like 
Aldebaran ( |van Leeuwen and Evans|1998|l, or e ven nearly ex 



tinct pu lsar tails (Roman ova et al.| (e.g. |2005] l). IBEX data 
( |McComas et al.|[2009[ |2011| i show that the Solar wind in- 
teracts with the ISM, and emits a steady stream of low en- 
ergy neutral atoms, producing a ribbon in the sky, completely 
unanticipated by theory and earlier observations. Here, we do 
not address these small scale features but focus on the expla- 
nation of the large-scale anisotropy. 

It has been argued that we sit in a local bubble, the walls of 
which are m ade up of various old SNR shells. However, a de- 
tailed study (Mebold et al.||1998 i suggests that this apparent 
bubble is not a coherent figure, but rather a motley assembly 
of filaments, clouds and shells, assembled into a bubble only 
in perception. These features represent the local history of 
star formation, HII regions, the effect of other stellar activity, 
and old violent supernovae. 

2.3. Balance of effects on cosmic ray anisotropy 

We conclude that old supernova remnants leave traces of their 
flow-field for much longer times, than traces of their individ- 
ual cosmic ray contribution. Therefore, the population of cos- 
mic ray particles in its scattering is more strongly influenced 
by these flow fields, than by the inhomogeneity of the sources 



in space and time, unless the source is unusually young. The 
anisotropy given by the flow fields is usually not dipole-like, 
since looking into different directions in the sky will see the 
thick surface of last effective scattering usually in different 
old supernova remnant flow-fields. 

It follows that there is an optimum particle energy to detect 
anisotropics, and at any larger energies we begin to smooth 
over several irregular flow patches (which clearly shows that 
the diffusion approximation is no longer adequate in describ- 
ing anisotropics), so that the anisotropics become smaller. 
Also, going to significant larger scales we begin to see dif- 
ferent flow patches, and so the directionality would become 
uncorrected at larger energies and thus larger scales as seen 
by inverting the scale energy connection: E ~ A^ fp ; this is 
derived by inverting the Kolmogorov spectrum with the ex- 
ponent 1/3. The IceCube data are consistent with both these 
effects, as the anisotropy is weaker at higher energies, and 
also uncorrelated in direction. Finally we note that of course 
the Sun may have a peculiar velocity with respect to its en- 
vironment: This has been s hown by VLB I ob servations to be 
a small effect, < 20km/s ( Reid et aT]|2009| . However, we 
need to emphasize that the uncertainties are large in such sim- 
ple arguments. Thus, we need to match particle energy and 
expected SNR scale: The old SNRs of about 30 pc match the 
mean free path for about 10 TeV. This should then correspond 
roughly to the maximum of any anisotropy; at larger spatial 
scales and thus larger energies the anisotropy is smeared out 
across several old SNRs; at lower energies we reach portions 
of an old SNR, not giving the full amplitude. Another im- 
portant aspect is that in this model there is no expectation of 
symmetry, since in one direction and the directly opposite di- 
rection the surface of last effective scattering will touch dif- 
ferent old supernova remnants. 

H owever, how does this co mpare to the elaborate calculations 
of |Blasi and Amafo] (2012a bl? These authors use the diffu- 
sion approximation, which of course begins to fail when the 
scale of the inhomogeneities becomes similar to the scale of 
the mean free path of scattering itself, a case which we argue 
we have here. So, using the diffusion approximation increases 
the apparent anisotropy when the scale of old supernova rem- 
nants is reached, while treating the scattering directly it be- 
comes obvious that then the "thickness" of the surface of last 
effective scattering itself begins to smear out anisotropics and 
so decreases it. 

Another aspect is that counting nearby sources of cosmic rays, 
when old supernova remnants not just produce a flow field, but 
of course also slowly disperse the cosmic rays as a function 
of time, produces a gradient in cosmic rays, which would also 
give an anisotropy: However, the data suggest that the aspect 
is a signifi cant deficit, but not a significant excess (Abbasi 
|et al.|2012) , a nd th is is more easily explained as a Compton- 
Getting effect (1935 ), which can have either sign, but does not 
have to be symmetric, as noted above. 

A key assumption in this model is the use of a Kolmogorov 
spectrum for the irregularities across the entire energy range 
of Galactic CRs, from 10 9 eV to 10 18 e V, a point on which we 
agree with Blasi and Amato ( 2012a|b i. Only this assumption 
allows the match in length scales. 

The amplitude at around 10 TeV would then be of order three 
times what is seen for the motion of the Earth around the Sun, 
so about 10~ 3 at most, and would correspond to angular scales 
of about 60 degrees. At larger energies the anisotropics begin 
to get smeared out. This is consistent with what is observed. 
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3. CONCLUSIONS 

In summary, we note as others have done before, that a vari- 
ety of effects may contribute to an observed anisotropy. The 
anisotropy suggested here, due to the Compton-Getting ef- 
fect in the magnetized flow field of old slowly disappearing 
SNRs would give rather large angle scale fluctuations, since 
the mean free path is valid in all spatial directions, and so an- 
gular scales, as noted above should be of order 60 degrees. 
The effect is a sum of uncorrected monopole components. 
That feature distinguishes this proposal from most other pos- 
sibilities, which can easily produce small scales. The key here 
was to assume that a Kolmogorov spectrum gives the length 
scale for a particle energy, which is quite different from other 
possible scattering descriptions. 

To the measure that the old SNRs argued here contribute to the 
main flow field influencing the surface of last effective scat- 
tering still have some of their own contribution to CRs, one 
might expect that contribution to be slightly flatter in its CR 
spectrum, at those locations in the sky, where the anisotropy 



is most visible. 

The effect suggested here predicts that at energies apart by 
more than an order of magnitude the anisotropics should be- 
come uncorrected, unless we sit in a very large-scale substan- 
tial flow-field. Averaging over several typical old SNRs scales 
the anisotropics should decrease with energy. Finally, to first 
order the anisotropics near about 10 TeV should have angular 
scales of order 60 degrees. 
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